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Abstract. We investigate the out-of-equilibrium dynamics of the one-dimensional 



■ quantum Ising model after a sudden quench of the transverse magnetic field. While 

for a translationally invariant system the statistical description of the asymptotic order 
' parameter correlations after the quench can be performed in terms of the generalized 

^ . Gibbs ensemble, we show that a breaking of translational invariance, e.g., by perturbing 

[ the boundary conditions, disrupts its validity. This effect, which of course vanishes in 

the thermodynamic limit, is shown to be very important in the presence of disorder. 



1. Introduction 



The experimental advances in the manipulation of ultracold atoms have paved the way 
^ ■ to accurately test the long-time dynamics of quantum many-body systems [1] . Starting 
■ from the observation of the collapse and revival of the superfiuid order parameter 
after a quench from a shallow to a deep optical lattice [2], a number of beautiful 
experiments probed the coherent relaxation dynamics in interacting quantum systems, 
with essentially perfect insulation from the external environment [21 [3l 01 [5l El [TJ [8]. 
Some typical features of adiabatic or quenched dynamics, that are by far harder to 
be identified in actual condensed matter systems, were observed neatly. Among the 
others, they include the phase dynamics emerging after splitting a one-dimensional 
Bose liquid [5] , or the absence of thermalization in out-of-equilibrium arrays of trapped 
one-dimensional Bose gases [1]. 

From a theoretical point of view, this experimental progress spurred renewed 
interest on the non-equihbrium dynamics of isolated many-body quantum systems, a 
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topic that was previously addressed mostly as an academic question (for a review, see 
Ref. |9]). In this context, an issue that received a great deal of attention has been the 
connection between quantum integrability and thermalization (or its absence) in the 
dynamics of a closed many-body systems [9]. The simplest protocol to study it in detail 
is the so-called quantum quench, i.e. a sudden variation of one of the parameters of the 
system's Hamiltonian, where the presence or absence of ergodicity should be detected 
by studying the state of the system long time after the quench has taken place. For 
non-integrable systems thermalization is supposedly occurring at the level of individual 
eigenstates [10] (see also [H]), and ergodic behavior is generally expected to emerge for 
few body observables [12]; actually this scenario seems too broad to be valid in any of 
such circumstances and it has been invalidated in some remarkable situations [T3 | [T ^ [T5]. 
On the other hand, for integrable systems the situation appears to be much less universal: 
non trivial constants of motion generally prevent the system from thermalizing in the 
usual sense of standard statistical mechanics. However, it was argued [TH [17] that it 
is still possible to describe the asymptotic state in terms of the generalization of the 
Gibbs ensemble (GGE), first proposed by Jaynes to take into account all constants of 
motion [TSj . 

The validity of the GGE for a variety of integrable models has been tested both 
numerically and analytically [HI HH [201 IHl [221 |23l |2l] . However, since the GGE does 
not take into account possible correlations between different constants of motion in the 
initial state, it was hinted in Ref. [25] that it should fail as soon as the breaking of 
global symmetries (like translational invariance) makes such correlations non negligible. 
The purpose of this paper is to show that this expectation is indeed correct: it is 
enough to slightly perturb the boundary conditions of an integrable model (or to 
make it non-homogeneous) to make the GGE fail. We will show this on the one- 
dimensional quantum Ising chain, the simplest non-trivial example of an integrable 
strongly interacting quantum system. This model has been the subject of various studies 
in the literature [26l[271[28l[29l[30l[ni[32l[33l[Tl[3l]. In particular it has been shown 
semi-numerically that for small quenches the asymptotic behavior of the order parameter 
two point correlation functions displays qualitatively a thermal behavior [301 131] . On 
the other hand, as recently shown by Calabrese, Essler and Fagotti [34], the correlation 
functions of the order parameter (as well as of the transverse magnetization) in the 
translationally invariant case can be quantitatively extracted from the GGE, finding 
agreement with the numerical analysis [301 [SI]- Here we will provide evidence that, 
if translational invariance is broken, the GGE fails to predict correctly the asymptotic 
state. We will first show this by studying a quench in the quantum Ising chain with open 
boundary conditions: for any finite system we find a significant discrepancy between the 
GGE and the exact analysis, which of course tends to zero in the thermodynamic limit. 
These discrepancies become of course very important in the presence of inhomogeneities 
not restricted to the boundaries, as we show by explicitly considering the presence of 
disorder both in the couplings and in the transverse field. The Ising model in the 
presence of completely random couplings and fields has been widely studied, see, e.g.. 
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Refs. [33 ESI EH EHl ESI SDl E]; here we are more interested in the regime in which the 
disordered part of the couphngs and of the fields is small, as compared to the uniform 
part. 

The paper is organized as follows. In Sec. [2] we set our notations, introduce the 
model under investigation and show the effects of the breaking of translational invariance 
on the structure of the eigenmodes of the system. In Sec. [3] we then frame the model 
in the quantum quench framework, and define the quantities we are going to discuss 
subsequently. Our results are then presented in Sec. HI Finally we draw our conclusions. 

2. Model 

We consider an integrable spin-1/2 quantum Ising chain in presence of random on-site 
disorder both in the couplings and in the transverse field. The Hamiltonian for a system 
of L sites is given by 

L L 

where a" (a = x,y,z) are spin-1/2 Pauli matrices for the i-the spin, J and F 
respectively denote the nearest-neighbor coupling and the transverse magnetic field 
strength, while rjj,hj G [—1,1] are dimensionless site-dependent quantities accounting 
for any disomogeneity in the model, and e sets the strength of the disorder. In the 
following we will consider both Open Boundary Conditions (OBC), implemented by 
setting 1 + e TjL = {or equivalently by imposing a^+i = 0), and Periodic Boundary 
Conditions (PBC), by assuming af+i = af. Hereafter we use units in which h = ks = 1, 
and set J = 1 as the energy scale of the system. 

The Ising chain can be solved by first performing a Jordan- Wigner transformation 
of the spin-1/2 particles into spinless fermions a~ = exp(— ivr J2i=i ^hus mapping 

Eq. ([1]) into the following quadratic Hamiltonian: 

-HsiT) = -J2il + e (cJsVi + cjc,+i + H.c.) - 2F ^^(1 + e h,)c]c, 

+ {-lf^(l + evL){c[cl + c[c, + H.c}j , (2) 

Np = Ylj=i ^j^i being the total number of c-fermions, whose parity operator (—1)^-^ is 
conserved, since it commutes with ?^£(F). In the general non-uniform case, the system is 
diagonalized by means of a Bogoliubov rotation [37] , introducing the fermionic operators 

L 

7/^ = 5^Kyi+^;/.c]). (3) 

The L-dimensional vectors = (ui^, ^2^,, . . . , ul^) and = {vi^, V2f„ . . . , VLp), for 
/i = 1, . . . , L can be always chosen to be real, and satisfy the coupled Bogoliubov-de 
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Gennes equations: 




where A, IB are L x L real matrices (A symmetric and B antisymmetric) whose nonzero 
elements are given by Ajj = — r(l +ehj), Ajj+i = Aj+ij = —(1 +erij)/2, Bj^+i = 
— Bj+ij = — (1 + £ ?7j)/2. If PBC are used, then the site j = L+ 1 is identified with j = 1 
and the following matrix elements are also present: A^ i = Ai^^ = (—1)^-^(1 + eriL)/2, 
Ml^i = —Mi^L = (—1)^^(1 + er]L)/2. In the disordered case e 7^ it is necessary to 
diagonalize the 2L x 2L eigenvalue problem in Eq. (jl]) numerically, in order to write the 
Hamiltonian 'He(r) in terms of its eigenmodes: 

ns{r) = 2j2^,(^^l%-\^ , (5) 

where > 0. The ground state of 'HeiT) is the Bogoliubov vacuum state \ipo) such that 
7fi\^o) = (Vyu = 1, . . . , L), with an energy Eq = - Yll=i ^z^- 

Before starting the discussion of our problem, let us comment on the form of the 
normal modes for the homogeneous system with PBC and OBC, displaying a crucial 
difference which turns out to be quite important in the quench dynamics. In the PBC 
case the index /i is the momentum index k = ivr/L, iSvr/L, .... The system can be 
equivalently solved by first applying a Fourier transform on the c-fermions and then 
diagonalizing two-by-two matrices, each one defined in a {k, —k) subspace, a procedure 
which is totally independent of F. On the contrary, in the OBC case due to the 
requirement of vanishing wavefunction at the boundaries, the solution involves standing 
waves of the form sin(A;j) where, however, the momenta k obey a quantization condition 
explicitly involving the transverse field parameter F ^2] : 

sinfc(L + l) ^ _1 

sinfcL F ■ ^ ^ 

As a result, the /c- vectors are not conserved if the transverse field is changed, for instance 
by a quench Fq — F (see below) and the Hilbert space can only be split in two subsectors, 
with odd and even number of fermions respectively, while /c-vectors are clearly mixed 
after changing F. A similar situation occurs as soon as e 7^ 0, thus showing that 
the invariant [k, —k) subspace structure for the translational invariant situation is very 
fragile and suddenly disappears in presence of any small imperfection: It is then unlikely 
to be the case in experimental implementations of such model. 



3. Quench dynamics and the Generalized Gibbs Ensemble 

The simplest non-equilibrium situation that can be established and controlled quite 
feasibly is a quantum quench in the magnetic field strength, that is a sudden variation 
of the value of F in Eq. ([1]). At zero temperature and in absence of any interaction 
with the environment, this consists in preparing the system in the ground state \iIjq) of 
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'H^i^o) corresponding to a given transverse field Tq. At time t = 0, the field is abruptly 
changed to some different value T ^ Tq. One then typically follows the unitary evolution 
under the modified Hamiltonian: 

|V;(t))=e-^^^(r)*|^o). (7) 

For closed integrable models, the generic absence of thermalization in terms of usual 
statistical ensembles is ascribed to the existence of non trivial integrals of motion. 
The way to take into account their effect is to define a Generalized Gibbs Ensemble 
(GGE) [16], which explicitly takes into account all such constants of motion through 
a set of Lagrange multipliers A^, 

PGGE = ^e-^.^-^^ (8) 

where Z = Tr(exp[- ^Jfi])- 

In the case of the quantum Ising model, Eq. ([T]), the number = = 7^% of 
quasi-particles constitutes a set of integrals of motion: every state can be identified in 
terms of the occupation numbers of all quasi-particle modes. The average occupation 
in the initial state can be computed by using the mapping from the fermions 7^°, which 
diagonalize ?ie(ro), to the fermions 7^ which diagonalize ^^(r): 

n, = iHll^llli^o) (9) 

L 

= E < ""^^ + € ""l """l" + <M < ^fu + <M <° 

jlu=l 

One can then construct a GGE starting from these conserved quantities. The value 
of each Lagrange multiplier is fixed by requiring (7^H^)gge = n^j,, such that 
= ^AjqiY" Hereafter we denote the expectation value of any observable O in the 

GGE with (a)GGE = Tr[pGGEO]. 

Notice now that for the quantum Ising chain the GGE takes the form of a 
generalized grand-canonical ensemble with a fi dependent chemical potential. In 
particular, all correlations between the occupations of states with different /i possibly 
present in the initial state are washed out in the construction. While in the ordered 
PBC case the {k, —k) structure of the solution is preserved by the quench, so that, for 
instance 

ii^olil^fM = a k'^k, 

with OBC and, more in general, when the system is non homogeneous, the mixing of 
the eigenmodes leads to 

(V^o|7i^7ri^o) 7^0 V/i,z/. 

As we shall see later, this reflects itself in a breakdown of the validity of the description 
of state ([7]) by means of the GGE. 
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3.1. Observables 

In the following we will concentrate on the asymptotic behavior of the two-point 
correlators of the order parameter o"^. We will examine both the equal-time correlator 
C^^(0, r) of two spins inside the chain at distance r, and the time-evolved autocorrelation 
function of the central spin, C^^(t,0). For the out-of-equilibrium system after the 
quench, they are explicitly given by 

C^Q^it^r) = {^{to)\ al/,^[H]it) mto)) , (10) 

where 0[H]{t) = QiHeir)t q ^-iHe{r)t (jg^otes the operator O in the Heisenberg 
representation at time t, while > is a lapse time after the quench, that is required 
for the excited quasiparticles to propagate along the chain 0. 

As shown below, for the ordered system (e = 0) the correlation functions (fTOj) decay 
exponentially both as a function of space and of time: CQ^{t ^ tQ,r) ~ e~^^^'^ and 
CQ^{t,0) ~ e~*''^<3, the correlation length and time Tq are determined, qualitatively 
and to a good accuracy even quantitatively, by the canonical ensemble average, that 
is only through an effective temperature set by the initial energy and the final 
Hamiltonian [30]. The effective temperature, Tes = (3~g is obtained by comparing 
the energy of the initial state with respect to the quenched Hamiltonian '^^(r) 
to the average energy of a fictitious thermal state in an effective canonical ensemble 
{i/jo\ne{roMo) = {HeiX)),an, where {6), an = Tr (Oe-z^^a^^ ) /Tr (e-'^'^s^-' ) denotes the 
canonical ensemble average [30]. Notice that while this effective description in terms of 
an effective canonical ensemble is able to reproduce the gross features of CQ^(t,r), its 
fine details require a more detailed analysis, which for the translational invariant case 
(ordered system with PBC) has been shown to be possible using a GGE [32l [34] • In order 
to assess the degree of validity of GGE in the presence of open boundary conditions, 
and eventually disorder, we will therefore compare the asymptotics of EqJTO] with 

CgGe(^, r) = (a2/2,[H](^) '^2/2+r)GGE • (11) 

Notice that the lapse time to has no effects on the GGE, since the time evolution 
operators e~^'^^^^^^° and e*^=(^)*o cancel out inside the trace. Since the quantities defined 
in Eqs. ( ITOl) -( |TT|) are in general complex, in what follows we will show results for their 
absolute value. 

The computation of both the correlators CQ^{t,r) and Cgge(^5''") can be easily 
implemented numerically, because the latter can be written in terms of the fermions 
(7/1) 7m) as the square root of a Pfaffian [26l [42], see also [37l |38l l40] . 

I The lapse time to is chosen to be sufficiently large such that CQ^{t,r) at given values of t,r has 
already reached a steady state, which is independent of to. Conversely, at fixed system size L, the lapse 
time to has not to be chosen too large, due to the appearance of revivals for very long times. 
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4. Results 

Different boundary conditions 

The simplest possibility to break the translational invariance (at least the most harmless, 
since it corresponds to adding one local defect in the transverse coupling between two 
spins, and it becomes irrelevant in the thermodynamic limit) is to slightly perturb the 
system by changing its boundary conditions. We follow this choice and discuss the 
changes in the asymptotic description of the order parameter correlator between the 
homogeneous PBC system and the OBC system. 

Let us begin with the spatial decay of the two-point correlators |C^^(0,r)|, and 
study the discrepancies between the quenched dynamics, |CQ^(0,r)|, and the GGE 
averages |Cqqe(0, r)|. As it is apparent from the upper panels of Fig. [1], the correlation 
functions always decay exponentially with the distance in both cases, and both for 
OBC and PBC: |C^^(0,r)| ~ e"*/^, where ^ denotes the correlation length. However, 
while for PBC there is substantial agreement between the asymptotic correlator and 
that predicted by the GGE (as it is apparent from the right panel), this agreement is 
definitely disrupted for OBC (left panel). In this case, contrary to the quenched case, 
the GGE data exhibit a dependence on the size L; this introduces some discrepancies 
with the out-of-equilibrium case (dashed line). The size-dependence is highlighted in 
the lower panel, where we plotted the correlation length and ^gge as a function 
of the inverse size \/L. Summarizing, the values of after a quench are essentially 
independent of the size and the boundary conditions (for the case presented here with 
Fq = 0.8 and F = 0.6, for example, we found ~ 49.5) and agree well with Cgge,pbc; 
vice- versa, ^gge, obc displays a systematic discrepancy which apparently vanishes in the 
thermodynamic limit L — oo. Ultimately, the discrepancies found between GGE and 
quenches with OBC are due to the fact that the initial state of the quench and the 
Hamiltonian after the quench are governed by different sets of A;-values, as remarked in 
Sec. [21 

We now turn to the on-site correlator |C^^(t,0)| as a function of time. As for the 
spatial decay, also in this case we observe an exponential decay: |C^^(t,0)| ~ e"*/"^*^, 
until a time t* where finite-size revivals appear[§|. This occurs independently of the choice 
of the boundary conditions, and both for the system after the quench and for the GGE. 
The differences in the GGE average between OBC and PBC emerge clearly also here, 
and closely resemble what we already observed for the spatial correlators. In the first 
case (upper left panel), the size dependence of |Cgge(^'0)I evident. On the contrary, 
the PBC-case is basically insensitive to L (upper right panel) [jj]. The out-of-equilibrium 

§ A rough estimate of the revival times is given by t* ~ N/vr, where vr is the maximum phase 
velocity of the spin chain over all the normal- mode velocities. Given the dispersion relation E^, then 
vr = Maxj,(9fc£'fc)|j.. For the uniform Ising chain with PBC, the maximum can be easily evaluated to 
give vr = 2Min(r, 1). 

II The computation of the autocorrelation function after a quench with PBC is not straightforward. 
The difficulty lays in the fact that '^j[H]i^)'^j connects states with different c-fermion parity which are 




Figure 1. (color online). Upper panels: expectation value of the equal-time correlator 
|C^^(0, r)\ as a function of the distance r from the central site, for OBC (left) and PBC 
(right), after a quench within the Ising ferromagnetic phase from Fq = 0.8 to F = 0.6. 
We show data both for the GGE averages and for the quenched dynamics. The colored 
lines denote GGE averages CQQ^{0,r) performed at different sizes L, as shown in the 
captions. The black dashed lines indicate Cg^(0, r) for L = 400 sites, and after a time 
lapse to = 100. Lower panel: comparison between the correlation length ^gge obtained 
with OBC (filled symbols; black circles for the GGE and red squares for the quench) 
and with PBC (empty symbols; blue triangles for the GGE and green diamonds for 
the quench), as a function of the inverse size 1/L. The dashed black line is a linear 
fit of the OBC data for the GGE at large sizes. Error bars have been estimated by 
changing the fitting interval over the data for r < L/2. 



correlation function CQ^(t,0) after a quench appears to be essentially independent of 
the size and of the boundary conditions, for times shorter than the revival time t*. If 
we compare the decay rate Tq (in the case highlighted in the bottom panel of Fig. Ej 
where we presented data for a quench from Fq = 1.6 to F = 1, we found Tq ^ 2.48) with 

therefore subject to different boundary conditions in the PBC case, see the last term in the r.h.s. of 
Eq. (I2])- We overcame this point by evaluating a c-parity conserving four-fermion correlator and then 
took its square root |26j . 




Figure 2. (color online). Upper panels: expectation values of the autocorrelation 
function |C^^(t,0)| on the central site as a function of time t, computed for different 
sizes L, with OBC (left) and PBC (right). All the curves denote GGE averages, except 
the dashed black line which is for the quenched dynamics. A quench from the Ising- 
paramagnetic phase Fq = 1.6 to the critical point F = 1 has been performed. Lower 
plot: comparison between the decay rate Tq and Tqqj, both for OBC and PBC, as 
functions of the inverse size 1/L. An estimate of the error bars is obtained similarly 

to Fig. m 



the rate Tqq^ extracted by averaging over the GGE, we found a quantitative agreement 
only for PBC, while some discrepancies are present for the decay rate averaged with 
open boundaries. In particular, at a given finite size, Tq q^^^ 7^ ''"gge obc sven if the 
agreement improves for increasing sizes. Despite the fact that we cannot be conclusive 
with only numerical results, for L -> 00 the results appear to have a tendency towards 
a convergence to a unique decay rate equal to the one for the quenched dynamics in the 
thermodynamic limit, where boundary conditions should be irrelevant. 
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Figure 3. (color online). Spatial decay of the equal-time correlator |C^^(0, r)| in 
the disordered Ising model. The various panels are for different values of the disorder 
strength e = (upper left), 0.025 (upper right), 0.1 (lower left), 0.5 (lower right). The 
dashed black lines are for the out-of-equilibrium system after a quench from Fq = 0.8 
to r = 0.2 (here we set a lapse time to = 300), while the continuous red curves denote 
the corresponding GGE averages. A system of i = 400 sites has been considered. We 
performed averages over 200 and 25 disorder realizations, respectively, for the quenched 
dynamics and for the GGE average. 



4-2. Effect of small disorder 

There are many other possible ways of breaking translational invariance. Disorder in 
the system destroys such invariance everywhere, in a much more substantial way than 
by changing the boundary conditions. Here we discuss what happens if a random 
component in the couplings and the magnetic fields is added, e ^ in Eq. ([1]). 

As it is apparent from Fig. |3l where we analyzed the spatial decay of the two-point 
correlator of the order parameter, the differences between the quenched dynamics and 
the GGE averages are striking. First of all we note that \CQ^{0,r) \ tends to a constant 
value for r — )■ oo, while IC'ggeIO' decays exponentially, irrespective of the value of 
e (remarkably, the effect of the noise in the GGE average is basically negligible, on 
the scale of the figures). This clearly shows that the validity of the GGE suddenly 
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ceases, as soon as a tiny amount of disorder is added in the system. In the panels of 
Fig. [3l we observe a quantitative agreement of the GGE average only in the first panel 
which corresponds to the limiting case of zero disorder. The differences between the two 
curves for r > 100 are due to finite-size effects. We remark that, in order to minimize 
fluctuations due to the disorder, we performed averages over different noise realizations. 
We found that, while for the quenched dynamics this was essential, the GGE averages 
are basically unaffected by the disorder: already with a single realization of the couplings 
a clean exponential decay clearly emerges (see the red curves in the figure). 



5. Conclusions 



We explicitly tested the predictions of the generalized Gibbs ensemble in the ID quantum 
Ising model. We have shown that the GGE is able to provide a quantitative description 
of the asymptotic behavior of the order parameter correlation function only in the case 
in which there are no correlations between the different constants of motion. This 
corresponds to the homogeneous, translationally invariant, situation pi]. On the other 
hand, as soon as the translational invariance is broken, the GGE fails to apply. 

In particular, we performed two different analyses. First, we added a local 
perturbation, by changing the boundary conditions from PBC to OBC: in this case, 
even though the GGE still exhibits a qualitative agreement with the out-of-equilibrium 
dynamics, a quantitative prediction is recovered only in the thermodynamic limit 
L -> oo. Second, we perturbed the system globally, by adding disorder both in the 
interaction coupling strengths and in the transverse field: in this case, the effects of 
the translational invariance breaking are much more evident: even a very weak disorder 
leads to a complete failure of the GGE in describing the quench dynamics. We found 
that, for the dynamics after the quench in presence of disorder, it is no longer possible to 
define a typical relaxation time, in accordance with what recently noticed in Ref. |13]. 
This is due to the localization effects, playing here a major role. Quite surprisingly, 
the GGE predictions are barely affected by such effects, thus exhibiting a qualitatively 
different (exponential) decay of the correlation functions. 
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